function x = precondU(b)

global U;
x = U \ (U' \ b);